FF No. 602(C) 


MSA 


X-641-66 tt %0k 

Revised November c .2 


u 


* r 


r. 


t w - h i i.a r 


5 


THE CALCULATED SPECTRUM OF iNVE f ' cc 
COMPTON SCATTERED PHOTONS B 


D 

fc-* 

Ln 

ru 

X- 

X- 

fc-» 


1 o 

: z 

! C 

s 

2 

30 

■< 

2 
a 


FRANK C. JONES 


n fi a - 1 1 s 5 3 

(ACCESSION DUMBER) 

^ / 

(PAGES) 



(THRU) 

(C^DE) 


(NASA CR OR TMX OR AD NUMBER) (CATEGORY) 



NOVEMBER 1967 


( 


LOAN COPY: RETURlf TO AFflh 
TECHNICAL LIBRARY* KIRTLAND APB, EJ 


GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 


REPRODUCED BY 

NATIONAL TECHNICAL 
INFORMATION SERVICE 

U. S. DEPARTMENT OF COMMERCE 
SPRINGFIEIO. VA. ZZ161 



TECH LIBRARY KAFB, NM 



X-641 -66-379 
Revised November 1967 


THE CALCULATED SPECTRUM OF INVERSE 
COMPTON SCATTERED PHOTONS 


F rank C. Jones 

Laboratory for Theoretical Studies 


November 1967 


GODDARD SPACE FLIGHT CENTER 


Greenbelt, Maryland 


PRECEDING PAGE BLANK NOT FILMED. 


THE CALCULATED SPECTRUM OF INVERSE 
COMPTON SCATTERED PHOTONS 

Frank C. Jones 

Laboratory for Theoretical Studies 
ABSTRACT 

We consider an electron of a given energy moving in a mono- 
energetic, isotropic radiation field. The energy spectrum of the pho- 
tons that are scattered by the electron has been calculated both exactly 
and in a greatly simplified approximate form suitable for astrophysical 
calculations. The approximation may be derived either by expanding 
the exact solution in a small parameter and keeping only the leading 
terms or by employing a simplifying physical approximation at the be- 
ginning of the calculation. The approximate spectrum is similar to one 
previously derived by Ginzburg and Syrovatskii, the principal difference 
being that the present one does not break down if E> (m e c 2 ) 2 where 
fta-'j is the initial photon energy and E the electron energy. We indicate 
the astrophysical applications of our approximate spectrum by calcu- 
lating the spectrum of photons scattered by electrons with an inverse - 
power-law energy distribution. 
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THE CALCULATED SPECTRUM OF INVERSE 
COMPTON SCATTERED PHOTONS 

I. INTRODUCTION 

In recent years the process referred to as inverse Compton scat- 
tering has had a revival of interest among astrophysicists. It was 
introduced in 1947 by Follin* as a mechanism for the loss of energy of 
cosmic ray electrons and was investigated by Feenberg and Primakoff^ 
and by Donahue^ in this context. Since that time it has been employed 

A Q 

in many treatments of cosmic ray electrons and the process was 
investigated in some detail by the present author in an earlier paper. 9 

It was first suggested as a source of energetic photons by Savedoff^ 
and by Felten and Morrison^ and has since received considerable at- 
tention^"^ from this point of view. Most of the calculations of photon 
spectra to date have been based on a rather simple approximation. It 
has been noted that the average energy transferred to a photon in a 
Compton collision is proportional to the initial energy of the photon and 
the square of the electron energy. This dependence on the square of the 
electron energy is reminiscent of the synchrotron process and for this 
reason the radiated photon spectra for a single electron energy is ap- 
proximated by a delta function spike at the average radiated energy. 

This spectrum is then folded into the distribution of electron energies 
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to produce the resultant photon spectrum* Although this method is 
known to give satisfactory results for synchrotron spectra it is now 
known 20 that the inverse Compton spectra are sufficiently different to 
raise some doubts as to its applicability in this area. However, it can 
be shown 21 that in the case of inverse power law distributions of elec- 
tron energies the method is applicable to both cases in spite of their 
differences. 

In the present paper we derive exact formulas for the scattered 
photon energy distribution for the case of an electron of energy y - E/mc 2 
moving through a region of space filled with a unit density of photons 
distributed isotropically with initial energy a x - Raj j / me 2 . We shall 
also derive several approximate formulas and discuss their validity 
in the light of the exact formulas. Similar calculations have been 
recently published by Baylis et al ,^ however, we find that our results 
disagree with theirs in several respects. In particular we disagree 
with their conclusion that a particular approximate spectrum is of as 
wide a validity as they claim. On the contrary we derive correction 
terms that become significant when certain conditions of validity first 
stated by Ginsburg and Syrovatskii^® are violated. 

In Section II we derive an approximate spectrum based on a simpli- 
fying physical assumption. The breakdown of this approximation will 
also be discussed from a physical point of view. In Section III the 
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scattered spectrum will be calculated exactly and compared (as well 
possible) with the approximate spectrum. We will see that the exact 
form is often not too useful for computation and the reason for this will 
be discussed. In Section IV we shall exhibit a method of expanding the 
exact formula in the small parameter that causes the trouble. This 
provides not only a method for computing the exact spectrum but pro- 
vides a systematic way of rederiving our approximate formula along 
with some correction terms. In Section V we will discuss some astro- 
physical implications of these results. 
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II. AN APPROXIMATE SPECTRUM 


In this section we shall derive an approximate spectrum by making 
a simplifying physical assumption concerning angles. Figure 1 illustrates 
the angles involved in a scattering problem as seen in the rest frame of 
the electron (E. R. frame). All quantities as measured in this frame 
will be primed and energies of the photon before and after collision 
(a.j and a respectively) as well as the electron energy y are understood 
to be in units of the electron rest energy me 2 . The polar angles Q ^ and 
9 ' are measured with respect to the electron velocity /3 = v/c (strictly 
with respect to minus /3; a l ’ /S ~ -o^ fi cos 8 ). The scattering takes 
place through a polar angle x' and an azimuthal angle <t>‘ where the /3 
plane is chosen as the 0' = 0 plane. 

The presence of so many angles along with the constraining rela- 
tions between them complicates the problem as we shall see later. It 
would greatly simplify things if we could eliminate some of them. To 
this end let us examine the angular distribution of the incoming photons 
in the E. R. frame. 

For photons, isotropic and monoenergetic with energy in the lab 
frame, the angular distribution in the E. R. frame is given by 


n' i) d^cos#^) 


d (cos 9 j'J 
2 y 2 (l - fi cos Sj') 2 


(1) 
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If yS ~ 1 half of the photons have polar angles within the range 0 <6^ < 9 ^, 2 
where 8 1 ' /2 % 1/y. In other words as the electron becomes more and 
more energetic the incoming photons appear to be more and more like 
a monodirectional beam with 6 ^ = 0. 

The approximation to be made is now obvious; we shall consider 
the electron to be energetic enough so that we may take 6^-0. A 
glance at Fig. 1 shows that in this case x ~ S' and since the scatter- 
ing cross section is independent of the azimuth 4>' we really have only 
one angle left to worry about. 

The energy distribution of the photons in the E. R. frame before 
collision is given by 

a i 

n' (dj') da^ = ^ 2 ^ ( a i* ’ a l/ 2 ’ y ’ a i ^ da i 


for 7 >> 1 and where S(x; a, b) is the characteristic function of the 
interval a , b i. e. 


S(x; a, b) - 1 for a<x£b 


0 for x < a , b < x 


The Klein-Nishina cross section for Compton scattering is given by 


Cr(a ' , dj' , y') - 


r ° 2(1+y,2) f, °- ,2(1 - y,)a i %l . ,i ■ , 

(3) 


2 [l +a 1 '(l-y')] i 
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where y' = cos X ' , r 0 = e 2 /mc 2 and f(a/, y') a/ [l + a/ (1 


The number of collisions per unit time t' is just N' ccr and since 


dN/dt = y ' 1 dN/dt' we have after integrating over 4>' 


7Tr o c 1 + y ' 2 


dtda/da'dy' 2a/ y 2 [l + a/ ( 1 - y ' )] 


a/ 2 ( 1 " y' ) 


(l + y' 2 ) [i +a i' (1 _ y')] 


:/ s(a' - f(a/, y ')] S(a/ ; a^Try, a x 2y) (4) 


Since da' da/ dy' = [l + a/ (1 -y ')] 2 da' dy' df we may integrate 
f immediately to obtain 


dt da' dy' 


nil. J (lty .« + » ,2(1 

2a,V Y 1 ' a ‘ a-y ) 


1 -a' (1- y') Vl -a' (1 -y 




We may relate a' to the final lab frame energy a by the Doppler 
shift formula a' = a/y(l -/3y' ). If we introduce the variable 17 = (1 “ ftf' ) 


we have 


77r 0 ca 


dtdadT] 2 y 4 a* ( 1 -a/y) 


T ) 1 ~ 27] + 2 + 


(a/y ) 2 1 gfo» 
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where we have assumed 1 - y' ^77 and where t ] 1 - a^2aj y 2 ( 1 - a/y) 
and 77 2 = 2 <z/a 1 (1 - a/y). The integral may be readily performed to give 


" _ o , _ 1 _ ( a/ y) 2 

V 21 nr) ^ ( 1 - a/y ) 77 

-*L 


(7) 


where the upper and lower limits, U and L depend on what part (if any) 
of the interval 77 x , r ] 2 lies within the limits 1/2 y 2 and 2. 

For a l ^ 4 y 2 < a <a t we have r ) 1 < 1 /2y 2 and 1 /2y 2 £ r ] 2 < 2. We then 
have, neglecting terms of order l/y 2 or less when compared to unity. 


d 2 N 
dt da 


7Tr o 2c Uy 2 a 

2y 4 a x \ a i 


(8) 


For a x £a £ 4 a } y 2 /(l + we have t 7 2 > 2 and l/2y 2 £T7 1 <2 and 


d 2 N 
dt da 


^ 7Tr 0 2 c 
a j 7 2 


2q" lnq" + (1 + 2q") (1 - q") 


1 .. 

+ 2 (1 +4a x yq") C 1 ” <1 > 


(9) 


where 


4a j y 2 ( 1 - a/y) 


and 


4y 2 


< q" < 1 
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In the above equations we see that the maximum value that a/y can 
have is (a/y) max = 4a j y (l + 4a j y) < 1 . If 4a t y << 1 then a/y << 1 and 
we have q" % a/4 y 2 . The last term in the square brackets may be 
dropped and we are left with the approximate spectrum of Ginzburg 
and Syrovatskii.^0 

Expression (9) is valid, however, no matter how large 4a t y may 
become. 

However, we must not assume that this approximation is uniformly 
valid. In fact it turns out that it is not a good approximation for 

a/a t as l/4y 2 . 

To see the reason for this let us consider what would happen if we 
took our assumption that 6 ' =0 seriously and transformed our photon 
energy distribution, Expression (2), back to the lab frame with no 
scattering at all. We would obtain a spectrum given by 

n(a) da = (a/a t 2 ) S(a; a 1 /4y 2 , a t | da (10) 

We can see from this that the approximation alone tends to populate the 
region of the spectrum from a^j4ry 2 to a x with no scattering at all and 
that this region will be exaggerated for small angle scattering as well. 

In other words, due to the large sensitivity of the Doppler shift formula 
to slight changes in 6' for small 9' neglect of these small deviations 
of 9‘ from zero introduces considerable error for small angle scattering. 
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Small angle scattering would be, of course, angles of the order of 
1 iy or smaller or for y' > 1 - l/ 2 y 2 . Since for most values of a/cxj there 
is a contribution from a considerable range of y' other than the region 
1 _> y' > 1 - l/2y 2 this error will be negligible. However, for the very 
bottom of the spectrum a/a 1 ^ l/4y 2 the contribution is entirely from 
the region 17 = 1 -/3y' ^ l/2y 2 and here we would expect a significant 
error. This will be borne out by the results of Section IY. 
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III. EXACT CALCULATION 


In this section we shall calculate in closed form the scattered 
photon spectrum for the case of an electron of energy 7 moving through 
a region of space filled with a unit density of isotropically distributed, 
monoenergetic photons of energy a^. As in the last section we will 
first transform the incident photon distribution to the E. R. frame and 
formulate the problem in this frame. The final energy will then be ex- 
pressed in its lab frame value a and we will integrate over all available 
angles holding a lf 7, and a fixed. 

In the E.R. frame the incident photon distribution has the form 


a ' 8 [a ' y(l -/3x' ) ~a x ] 

n' (a ' x' da,' dn' = ; — da.' dfJ' 

' 1 ' 1 a i 4 tt7(1 -fix ' ) 1 

where x' = cos#' and dfi' is the element of solid angle 2 ? 7 dx'. Expres- 
sion ( 11 ) may be obtained by noting that x) da 1 dQ is a density and 

hence transforms like an energy. If we divide by the energy a l we ob- 
tain an invariant and hence 


n (aj , x) da 1 dfi - ~ n' (a^ , x') da^ dfi' 


and 


n' x # ) da i dfl' 



n (a x , 


x) da x dfl 
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Expressing a x and x in terms of and x' completes the derivation. 




The cross section is given by expression (3) 


o-jaj', a', y') dfi' (y')da' 


r„ 2 (1 + y' 2 ) r a i' ! (i -y' ) 2 

1 + 


2 i 1 +a i' (i - y')] \ ( J + y' 2 ) [i +a i (i- y')]J 


S (a' - f (a x ' , y')j 64>‘ dy ' da' . (3) 


Since dN/dt - n' ca-/y we have 


d®_N 

dt dx' dy' dc f>' daj' da' 


r o 2 c s ( a i r(i-/3x) ” a i) 
4y 3 ( 1 - fix' ) 2 


5(a' - f (a j' , y)) 


X 


i +y' 2 

a i' 2 (i - y') 2 

r , . ,2 

r # t 3 

[l + a j (1 ~ y )] 

[l+a; ( 1 ~ y )] J 


(12) 


where we have used the relation - -y(l ~/3x'). Once again em- 

ploying the relation da x ' = [i W (l-y')] 2 df we may integrate over f 
immediately. We also note that since z = cos 9 9 - cos 9 ^ cos \ 9 
+ sin sin x* cos 0' we have 


d0' 


2dzJ 

(l_x'2 -y'2 _ z »2 +2x / y - z * J 1/2 
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With this substitution we now have 


d 5 N 

dt dx' dy' dz ' da' 


cr, 


2 


a' y(l - fix’) _ 
_1 -a'(l ~y') 
y 3 ( 1 - /3x ' ) 2 


a 


1 


x 





(13) 


Transforming a' to a using the relation 


a 

7(1 - /3z 1 ) 


we have 


d 5 N 

dt da dx ' dy 9 dz 1 


cr 0 2 -y(i-y') 

2 y 4 (1 -fix' ) 2 (1 - /3z' ) 


c 


S a( 1 _ /fix' ) 


a i (1-/8*') 


a. a 

+ — (l-y')| 


l +y 


< 2 


a 2 (1 - y' ) : 


• 1/2 


7 


( 1 - /3z ' ) l-/3z' -y (l-y*) 


(14) 


where J = 1 - (x' ) 2 - (y' ) 2 - (z' ) 2 + 2x' y' z' and where we have made 
use of 


S(Ax-B) 
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From this point on the object is to integrate over all possible 
values of x', y', and z' holding a, y and fixed. For a given set of 
values for the parameters a, y and only a certain volume of x' , y', 
z' space (possibly zero) will be compatible kinematically with this 
particular choice. This requirement is expressed by the condition that 
the Jacobian of the transformation from $>' to z' be real or that J > 0. 
Inspection of the form of J shows that the requirements that |x* | , ly' I 
and | z ' | be 1 1 is automatically fulfilled by keeping J>0 unless all 
three variables are simultaneously out of bounds in such a way that 
(x' y' z' ) >0. Therefore this requirement need be consciously en- 
forced on only one of the three variables. 

It is immaterial which order we choose in integrating the three 
variables and we arbitrarily choose the order x', y' , z'. The first 
integration is trivial because of the delta function and we obtain 


d 4 N 

dt da dy ' dz ' 


a 

y 4 a 2 

l-/3z'-f (1-y') 

( 1 " fiz ' ) 


where 


x (1 +y' 2 + 


(a/y ) 2 (1 ~ y' ) 2 


( 1 - /3z ' ) l-/3z'““(l-y') 


(/3 2 J)- 1/2 (15) 


/3 2 J = (/? 2 + e 2 + 2/3ez') (y 2 -y') (y' -yj 


(16) 
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(17) 


y 2 - y 0 + 8 
y x = y 0 ■ 8 


yo 


(V +pz' ) (p + ep- 1 + /3z* ) 
p(/3 2 + e 2 + 2/3ez ') 


(18) 


S 


yS(l - z' 2 ) 1 2 [p 2 /3 2 + 2pe(l - p) (1 - /3z') - (p - 1 + /3z' ) 2 ]^ 
_ ^(7 2 -£ 2 ’-2P6z') 


(19) 


and p - a/cij, e - a^/y. 

The integration over y* may be facilitated by the transformation 
y' = y 0 + St? where -1 <77 < 1. We then have 


d 4 N 


cr 0 «• 


/ 


y 0 2 + 2y 0 St ? 


+ 8 2 7] 2 


dt da dz ' dr] 


2y 4 a 2 


(/3 2 + e 2 + 2/3ez') 


'M 2 


V 


(1 - /3z' ) 


(a + bp) 


(q/y) (1 ~ y 0 - St?) \ 
(a + bp) 2 ] 


5 \ 1/5 


( 20 ) 


where a - l~/3z' (a/y)^l y 0 ) ; b ao/y. 

Integration of p from -1 to 1 gives 

d 3 N = 77r o 2 C(X / 1 

dt da dz ' ” 2y 4 a^ 2 (/3 2 + e 2 + 2/3ez ') 1/2 ((l-/3z') 2 


y 0 


(1 -0 z') (a 2 -b 2 ) 


1/2 


2y 0 y 2 y 0 ya 


ay 21 


a 2 y 2 


a(l-/3z') a(l-/3z')(a 2 -b 2 ) 1/2 a 2 (l-/3z') a 2 ( 1 -/3z' ) (a 2 -b 2 ) 


2\l/2 


a ( 1 _ y 0 ) 


a 2 S 2 


y(l -/3z')(a 2 -b 2 ) 3/2 y 2 (1 -/3z')(a 2 ~b 2 ) 


3/2 


( 21 ) 
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After some manipulation we have 


[( 1 - a/y ) 2 - 1/y 2 ] (l - /3z ') 2 + [2a/y 3 ] ( 1 - /3z ' ) 
(/3 2 + e 2 + 2/3ez') 


[y 2 (1 - ep) 2 ~ l] (1 -/fc* ) 2 +26/0(1 -/3z') 
y 2 (/3 2 + e 2 + 2/Bez') 


(22) 


If we introduce the variable £ - l-/3z' and the following quantities 
Ej = (/? 2 +e 2 + 2e)-2e£ = ( 1 + e ) 2 - 1/y 2 - 2e£ 


E 2 - y 2 [(1 - ep) 2 - 1/y 2 ] i 2 + 2epl, 


the final integration over z' may be performed in a straightforward 
manner, and after some rearrangement of terms we obtain 


d 2 N 
dt da 


ttt 2 ca 

2y 4 J3a 2 


[*•(£*) -r(U] 


(23) 


where £ ± are the upper and lower limits of the integration in £ and the 
function F is given by 


F(0 = f 1 (0 + f 2 (0 + f 3 (0 + f 4 (0 
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f,(C) = e; 1 -' 2 




1 + aa 


OLj - a 


a 


\ 1 ” oi/y) 2 - 1 /y 2 ] 


y 


+ i + T-U-r 


(25) 


MO = 


27' 


(y + a i ) 2 + aa i 


a 


(y + aj 2 - i 1 2 \ ( y +a i) 


cosh 


- 1 


(y + a.j y - 1 


2a ! 


1/2 


(26) 


-2y 2 


0 (O 


[o-a) 2 - 1] 1 


/ 2 


(y - a) 2 + acO 

1 + ~ 7^ H sinh 1 

(7 “ a) 2 - 1 


[(7- a) 2 - l] yC 


1/2 


~ 2 y ^ 


1 


a [l - (y - a ) 2 ] 12 l l-(y-a) 


(y - a ) 2 + aa ^ 

- 


sm 


-l 


2 a 

for y - a > 1 
1/2 


[l - (7 - a) 2 ] yi 


(27) 


2 a 

for 7 - a < 1 

We now turn to the question of determining the limits of integration 
£ ± * These are determined by the requirements that the quantity S be 
real and in addition that |z' j <1 for we see upon inspection of Expres- 
sion (19) that S may be real for certain values of z' that violate this 
condition. These two requirements are fulfilled if the quantity £ lies 
between the values 1 ± /3 called the boundary lines and simultaneously 
lies between the values 


i±(P) ~ p[(l+e-ep) ± j/(l + e - ep ) 2 - l/y 2 J (28) 

called the boundary curve. It is easy to see that at p = 1 the boundary 
curve intersects the boundary lines. At p s = l + (y-l)/ 7 € = 1 + (y - l)/a l 
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the radical in (28) vanishes and the boundary curve becomes imaginary. 
This clearly represents an absolute upper limit on p (the other real 
branch of the boundary curve for even larger p can be shown to lie en- 
tirely in a region of z' > 1). The physical significance of this limit is 
quite simply seen if we write it as a. = + (y — 1) . At this limit the 

scattered photon has picked up all of the electrons kinetic energy in 
the collision. 

This limit is not usually reached in any situation that will interest 
us since it only occurs when the initial photon momentum is of the order 
of or greater than that of the electron. Figures 2a, 2b and 2c illustrate 
the three different situations that can exist. It is quite obvious that the 
usual situation in astrophysics will be that depicted in Fig. 2a where 
for tXj _< + y( 1 + /3)J the maximum value of p is given by the point 

where the lower boundary curve intersects the upper boundary line; 

(p) = 1 +/3 or 

- Ill . 

Pc l-/3 + 2a j/ y 

In the relativistic limit 1-/3^ l/2y 2 and we have 

4y 2 4a i^ 2 

Pc ~ 1 + 4a x y or a max % 1 + 4aj y 
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which is just the result derived in Section II. The minimum value of p 
is always given by 

1-/3 

Pm l+/3+2a 1 /y 

which in the relativistic limit is 

p ^ 1/4 y 2 or a . ** a, /4y 2 

also a result of Section II. 

The formula given in expressions (23-27) are not very useful in 
most astrophysical applications either for insight, since they are quite 
complex, or for direct computation since they require that terms of the 
order of (y/a ) 2 be balanced out to yield a true leading term of order 
of y 2 dj/a. This requires computation to be carried out to an accuracy 
of (aaj p) % to obtain an answer that is correct to an accuracy of p %. 
Since the quantity a.a x can often be quite small direct application of 
Expressions (23-27) is in general not very satisfactory. 

In the next section we shall discuss various expansions of the func- 
tion F(£) which will be useful not only for computing the spectrum to 
any desired order but also for recovering a simple approximation with 
a wide range of validity. 
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IV. EXPANSIONS OF EXACT FORMULA 

The chief difficulty in computing directly with our exact formula 
1 - is the fact that the quantity ep = a/y is often quite small and appears 

as (ep)' 1 , (ep)" 2 in some of the terms. This in itself suggests the way 
out, namely an expansion in this or some other small quantity. The 
! quantity that turns out to be most useful as an expansion quantity is 

e = a^/y. This quantity is small in almost all physical applications 
and becomes smaller the more energetic the electron. 

The first step in this procedure is to expand the functions E ^ 172 , 
E 2 " 1/2 , cosh -1 , and sinh " 1 as power series in the quantities 2e(/[(l + e ) 2 - I/ 7 2 ] 
and 2ep/[y 2 (1 - ep) 2 - l] ( . 

The expansion in the first quantity can be easily shown to be con- 
vergent for all allowed values of the parameters and the second ex- 
pansion is convergent so long as we have the condition 

2j gTg _ 2 + 3/3 

a i < 7 + 9/3 " 7 ( 1 +/3) (7 +9/3) 

~ fZ/S - 5/32y for 7 » 1 

This will be true in most cases of interest and in those situations 
where it is not true the expansion will still converge as long as 

ep < 2 + /3 - 2 /Ti~3 ~ 0.172 - 0.146/7 2 , 7 » 1 
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If this condition is violated we see that ep is not a small number and 
there is no real need for the expansion. In the following we shall always 
assume that is small enough that this first expansion is fairly rapidly 
convergent. 

At this point we have an expression for the function F(£) of the 
form 


F(Q 



(29) 


where the term £°/0 = In £ arises as the leading term in the expansions 
of cosh -1 and sinh -1 . There are no terms independent of £ since all we 
are interested in is the quantity F (£ + ) -F(£_) and such terms would 
make no contribution. The coefficient K n (e) is a rather complicated 
function of e whose dominant term is of the order of e U" 1 !"2 > 

K n (e) contains the two denominators which are functions of e 

Dj = (1+e) 2 - 1/y 2 

D, = (1 - ep) 2 - 1/y 2 (30) 


These denominators appear in half powers of various orders (typically 
I n | / 2 ) and are the next items on the list to be expanded in e. Before 
proceeding, however, we must first decide how we are to order the 
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quantity p = a/a,. Recalling the limits on p, l/4y 2 < p < 4y 2 /(l +4a, y ) 
we see that the question hinges on the magnitude of a, 7 . If a, y « 1 
we have l/4y 2 <p£ 4y 2 and e « l/y 2 . We may then consider p to be of 
0(1) and expand in e and ep as well. On the other hand if a, y > 1 then 
e > l/y 2 and 0 ( e ) _< p < 0 ( l/e ) and the quantity ep ranges from 0 (e 2 ) to 
0(1). In this case we must expand in quite a different manner. 

We shall now consider the case where a, y << 1 and treat ep as 0(e). 
First noting that 

D, = ( 1 + e ) 2 - l/y 2 = (1 + l/y + e) (1 - l/y + e) 

we may expand the denominator as 

00 

D .'" 1 E p -.- ( I/y2 > (31 

n = 0 V 


where 


P n, m _ a m (P ) a m (q) (! _ i/y ) 15 (! + 1 /> , ) q 


(32) 


p + q-n 


and the a (p) are the expansion coefficients of 


1 ~ x ’ )_m _ a,n ^ xf 
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given by 


a m (p) - m(m + l)(m+2) * * * (m+p-1 )/p! (33) 

where a m (0) = 1 and a m (~n) =0. We may also write P n (1 /y 2 ) as 

/ 

P 

p n, m ( l/ y 2 ) " 2Z a,n (p)a2m+2p (n_2p) ^ 1/7 ^ P (32 ’> 

p = 0 

where p' is the largest integer not greater than n/2. 

A completely equivalent expansion exists for D 2 so our coefficient 
K n (e) may be in turn written as a power series in e of the form 



F m (£) is a rather complex function of £ but for completeness we shall 
give its general form for arbitrary M. First, however, we note that 
the expansion of F(£ + ) -F(£_) is obtained by simply replacing £ n by 
£+ n ~ £_ n wherever it appears in F„ (£)• what follows we shall use 
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the notation 


Z(n) = {i+~i-)/ n for r» t 0 

Z(0) - ln(£ + /U • 


We may now write the general expression for F M (£ + ) - F M (£_) _ F, 


M' 


/ . . V" 1 „ . „ I^M-N.N + 5/2 (l/7 2 ) T(N + 1) 

F m = (-1//S 2 ) M ^ (-2) N + 2 a 1/2 (N + 2)^ p- K — 1^2—^ Z(N + 3) 

n = -2 


+ - N 4 - 2 - Z(N + 1) - < N + 2) 2 (1 + Z(N + 2) + (-p) M 


P 


(> vj fp] 


M-N.N + 3/2 ( : l/y2 ) [2(N + 1) (N+ 2) (N + 1 + Np) (N + 2) 




(2N + 3) p 


Z(N) - 


( 2N + 3) p 2 


Z(N+ 1) 


+ ( _ ^) M (2N+3) ( 7 2) 


(N + 3 + (N + 1) y 2 ) Z (-N - 1)1 


N + 1 


^M-N.N + l/2 (^■ // ’> /2 ) 


yS 


(N + 1) (N+ 2) (1 + N/P) 

( 2N + 1) (2N+ 3) Z(N - 1) 


. N (N + 1) (N+ 2) /Z(l-N) , Z(-N) \ 

(~ P ) (2N + 3) ^ y 2 j n-x + ^ (y 2 ) N / 


+ (-p) M N(N + 1) (N + 2) y 2 


+ ^^M-N,N-l/2 ) 


(N-l)(N)(N + l)(N+2) Z(l-N) 

(2N-l)(2N + 3) V ) ( y 2)N-i 


(~P) 


( 35 ) 
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In deriving the above expression much use has been made of the recursion 
relations for the various a o (p) coefficients, i.e., 


1/2 


(N) = 


2N + 2 

(2N + 1) (2N + 3) 


l 3/2 ( N + !)* etc - 


Explicit calculation of the M - -2, -1 terms shows that they are identi- 
cally zero so the leading term in our expansion is of zero order in e . 
Our series may then be written 


F(C) -f(£.) = 2Z F “ e “ ■ 

M= 0 

It should be noted that there are terms with (-p) M appearing in F M and 
in the case p % 0(l/e) each term in the series will be 0(1) and our 
expansion breaks down completely. 

Although the expression for F M looks rather formidable, it is a 
straightforward matter to program a computer to evaluate our func- 
tion to any order in e that is desired. 

It is interesting to examine the zero order term F 0 . 

F 0 = f3~ 5 (2 + 2/(y/3) 2 + 3/2(y/3) 4 ) Z(-l) - /T 7 (2 + 3/(y^) 2 ) (l + 1/p) Z(0) 

+ /3" 9 (l+6/p + l/p 2 + ^J^ 2 ^Z(l)-/3- 9 (3(l + p)/p 2 ) Z(2) 

+ / 3" 9 (3/ 2p 2 ) Z(3) (36) 
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If we also assume 0» e we have £ + - 1 + 0 and to zero order in e , 


l 


p >_ 1 so 

that 

Z(-l) 

= {l+0-p(l~0j) 

Z(0) 

ln (y 2 (i + /6) 2 } 

Z(l) 

= l+/3-p(l-/3) 

Z(2) 

= \ (1 +/3) 2 - J P 2 (I"/?) 2 

Z( 3) 

= 3 C 1 + £> 3 P 3 (1'0) 3 


Expressions (36), (37) and (23) may be combined to obtain an approxi- 
mate spectrum that is valid for 0» e and <Xj y« 1. 

We may further simplify this approximation by assuming that the 
electron is relativistic, 0 * 1 , 1 -0 * 1/2 y 2 . We then have, neglecting 
terms of order 1/y 2 as compared to 1. 


d 2 N 
dt da 


27tt^c 
dj y 2 


{2qlnq + 0 + 2 q) ( 1 - q) + 0 (l/y 2 )} 


where q = a/ (4a y 2 ) . This is just the approximate spectrum of Ginzburg 

and Syrovatskii.20 


For p < 1, £ + = p(l + /3) and - 1 -/3. This gives 


Z(-l) = ~ [p(l+/8) - (1-/3)] 

Z(0) - In [py 2 (1 + /3) 2 ] 

Z(l) - p(l +/3) - (1 ~/3) 

Z(2) = l [, o 2 (1 + /3) 2 - (1-/3) 2 ] 

Z(3) = | [p 3 (1+/3) 3 - (1-/3) 3 ] (39) 

Once again neglecting terms of order 1/y 2 we have 


d 2 N 
dt da 



|(q' ~ 1) (1 + 2/q' ) - 2 In q' + 0 (l/y 2 )} 


(40) 


where q' = 4y 2 p - 4 y 2 a/a 1 . We note that this is just our approximate 
spectrum, Expression (8) with additional correction terms that become 
important at the bottom of the spectrum where a/a l ^ l/4y 2 . These 
correction terms are expected on the basis of the discussion at the 
end of Section II. 

We now consider the case where y is of order unity or greater 
and ep may become of order unity. We first note that if a 1 y>l then 
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1/y 2 < e and it would be inconsistent not to expand in 1/y 2 as well as 
in e . 

The expansion procedure is very much the same as before. The 
denominators D x and D 2 (Expression (30)) may be first expanded in 1/y 2 
and then the term (1 +e ) 2 that comes from D 1 is expanded in e . The 
term (1 - ep ) 2 arising from D 2 , however, is not expanded and all factors 
of p are combined with an e to make terms of 0(1). The resulting ex- 
pressions are then grouped according to the power of e and the power 
of 1/y 2 to give a double power series expansion of F(£) as 

00 

F(q - F(£.) = £ F n P (l/y 2 ) N e p 

N ,P- 0 

where F N p is given by 


N,P 


P 

y \ a 1/2 (P"p) a P ~p + 1/2 a 2(N+P-p) + l P 


P = 0 


x <Ml - 


(P-p)(P-p-l) (P"P + N) 

2(P +N ~ p/2 + 1/2) (P + N- p/2) (ep) 


.) 


Z(P-p- 1) 


<(P-p+2N+2) p _ p 


ep 


( e /^) 2 


Z(P-P) 


_ _ _ i (2N+2P-P + 1) (2N + 2P-p+2)\ 

+ '-^ 2(N + P ~ ~ p + ' i ~yTp j Z(P-P + 1): 


(41 
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1 


N 

_ V" 


a 3/2 (N -m) 


m - 0 


x 


PO 


l ( e /°) 2 ( i - e P ) 


Z(m - N + 1) 


+ 


[ S po [2(N-m) 2 + N-m-2] 

(ep) (1 ~ep) 


S P1 (1 + 2N)\ 


+ 



1 + 2N 

(1 - ep) 2 


( 1 + 2N - 2m) (N-m 


(1 - ep) 


)' 


+ S 


PI 


(1 + 2N) (1 + ep) 

(1 +2N - 2m) (1 - epy 


|Z(m-N- 1) 


(4 1 ) con’ t 


where 


S 1 if m - n 

m, n 


- 0 if m t n 


Once again we examine the lowest order term F 0 0 . 


2 ( 1 + ep) 

F o.o = Z(-l) + -Z(0) -y^y-Z(l) 

+ 2 . Z(0) + Z( :. ^ 

(ep) 2 (l-ep) ep(l-ep) (1-ep) 2 


_ Z(l) 
1-ep 


2Z(0) + 

(1-ep) 


( 


(±py 


P e P) (l-ep) : 


z(-i) 


(42) 
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To lowest order in e and 1/y 2 the boundaries of £ are given by C + 
= l+ /3^2;£_ = p/2y 2 (1 - ep) so we have 


'o,o (1 - ep) 


In 


i4y 2 (1 - ep), 


H 


4y 2 (1 - ep) 



(1 - ep)' 


1 ( 6 P ) 2 Ay 2 (1 - ep) 

2 (1 - ep) 2 \ P 


■) 


1 (43) 


Combining this with Expression (2 3) we have 


d 2 N 
dt da 


2' 7Tr Q c 


a i 7 


2q" In q" + ( 1 + 2q" ) ( 1 - q" ) 


1 ( 4a i yy") 2 

2 (l + 4a } yq") 


d-q") 


(44) 


where now 


4T 2 (1-ep) 4a 1 y 2 (l-a/y) 

We see that we have recovered our approximate spectrum of Section II 
Expression (9). Now, however, it appears in a complete mathematical 
setting as the lowest order term in a double expansion in € and 1/y 2 
where a/a x - 0(l/e). 
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We have seen that the assumption that a/a 1 » 1 is necessary in 
deriving this formula. We now may ask what happens Lkj this approx- 
imation when we no longer have a/a l = 0(l/e) but return to the region 
where ep - 0(e). We can see at once from Expression (41) that due to 
terms containing various powers of ep and denominators (1 - ep) to 
various powers that a term that was originally of a given order will now 
contain contributions of all higher orders in e. For any approximation 
to a given order this does not cause any loss of accuracy. What does 
hurt is the presence of terms containing (ep) -1 and (ep) -2 . This means 
that any given order now has contributions from terms that were 
previously as much as two orders higher. To maintain our zero order 
approximation we must now include those parts of eF 01 and e 2 F 02 that 
contribute to zero order in e. These terms may be found in a straight- 
forward manner and we find that eF 01 to zero order in e gives (l/b)(6Z( 1) 
- 3Z(2)- 2Z(0)) and e 2 F 02 gives 

\ (Z(1)-3Z(2) +f Z(3)J . 

These terms are the same terms in 1/p and 1/p 2 that appeared in 
Expression (36) (neglecting 1/y 2 ). They did not appear in Expression 
(38), however, since they are always at least (1/y 2 ) smaller than the 
leading term inZ(-l) which is 2y 2 /p. For that reason they should not 
be included in our present formula. 
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We, therefore, offer Expression (44) as an approximation to the 
spectrum of inverse Compton scattered photons that is accurate to 
zero order in e and 1/y 2 for values of a such that a 1 < a < 40^ y 2 /(l + 4a t y) 
and for a^as large as desired. Expression (44) is not an entirely 
consistent formula in that it always contains contributions from higher 
(and hence negligible) orders but it is complete in that it always in- 
cludes jill_zero order contributions. 

In the situation where p = 0(l/e) the exact formula should be ex- 
panded once again, this time considering e/p = 0(1). However, if we 
are interested only in the lowest order approximation a simple in- 
spection of Expression (35) will suffice. Keeping in mind that for p < 1, 

[(V p)"-l] 

z<N> = N(vr 

we see that for every value of M there are terms of order p’ 1 and 
higher but of no lower order. Therefore, Expression (35) to zero 
order in e and lowest order in 1/y 2 will give us our dominant term. 

This is exactly what we obtained in Expression (40) so we see that 
this formula gives the correct approximation to lowest orders in e and 
1/y 2 no matter what the magnitude of a^y. This could have been ex- 
pected from our discussion in Section II. 
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In Figure 3 we compare the approximate spectrum of Expres- 
sions (40) and (44) with a computer calculation of Expression (35) 
correct to order e 5 . We see that the electron does not have to be ex- 
tremely relativistic for the approximate spectrum to give a good 
representation. 
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V. ASTROPHYSICAL APPLICATIONS 


In astrophysics, inverse Compton scattering provides a mechanism 
for energy loss of high energy cosmic ray electrons and a source of x 
and y radiation whenever energetic electrons and soft photons exist 
together in a region of space. The energy loss effect has been calcu- 
lated exactly by the author in a previous publication.^ However, it 
would be of interest to see how well our spectrum, Expression (44), 
serves in giving the correct energy-loss formula. In keeping with the 
spirit of our approximation, we shall assume that Expression (44) is 
valid for 0£q"£l even though we know it is quite invalid for q" < l/4y 2 . 
When we consider effects that depend on the entire spectrum the region 
0<q"<l/4y 2 contributes a part that is 0(l/y 2 ) and hence, negligible. 

The energy loss is given by 



where c T - (8/3)77r 0 2 and b - 4a 1 % The indicated integrals are performed 
in a straightforward manner to give, after some rearranging of terms. 
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3 (J T c y 

b 2 


<j(b/2 + 6 + 6/b) In ( 1 + b) 
v. 


b 3 + 6b 2 + 9b 


(1 +b)" 2 


- 2 


+ 2Li 2 (-b) 


(46) 


where Li 2 is the dilogarithm. If we make the substitution b = 2a we 
have 



F(a) = 7 


(a + 6+ 3 /a) ln(l +2a) 


a 3 + 24a 2 + 18a +4 j ( 1 + 2a) 


-2 


- 2 - 2Li 2 (-2a)J (47) 

This expression may be directly compared to Expressions (13) and (14) 
of Reference (9). It can be seen that the present result is equal to that 
previously calculated (for the mono-energetic background case) if one 
sets e = 1/y 2 = 0 in the exact formula. 

In considering inverse Compton scattering as a source of x and y 
rays we are interested in the radiation from electrons with a wide 
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distribution of energies. In astrophysics the inverse power law is one 
of the most commonly occurring distributions so we shall consider 
the spectrum 



(48) 


dy 



dcf 

( pq ” ) 3/2 (1 + pq ") 1/2 


where p - aa l we may write Expression (48) as 


2 r 

r O) _ a (r+ 1 )/ 2 F (P’ r ) 


(49) 


where 


F(P. O 


i 


i q «(r-i)/2 


[2q lnq" + ( 1 + 2q" ) ( 1 - q" ) + 4pq" ( 1 - q" )] 


1 + 


00 " 


dq" (50) 
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It is easy to see that for p - aa t << 1, F(p, T) is essentially independent 


of p and we obtain 


R(a) 




+ 1 


if + 3 T r + 1 


a) 4' 


_ ^ jr-iv 2-(r+i )/ 2 

+ 5 (r + 3)2/ a - a 


= c' a-( r+1 > /2 (51) 

This is just the well known approximate spectrum of Felten and 
Morrison which they obtained by approximating d 2 N/dt da by a delta 
function S (a - 4/3 a l y 2 ) . It can be shown^ that this spectrum is a good 
approximation for a much more general assumption about the form of 
d 2 N/dt da than that made by Felten and Morrison. 

For p >> 1, on the other hand there is no expansion in powers of p 
or inverse powers of p that will be valid for the entire range of the 
integration in ql' However, the integral may be broken up into two 
pieces as 


J -i/p pi 

+ I 

0 /p 

and appropriate expansions made in each range. The leading terms in 
p may be extracted and after a certain amount of re -summing of 
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coefficients we obtain the assymptotic form of F(p, F) for p» 1 as, 


F(p, F) 


In p - c(T) 

2 r p( r+ i)/2 


( 52 ) 


where unfortunately 


c(r> 



(53) 


Inserting (52) in (49) we have for the case a >> 

R(a) ~ (in aaj +c(F)j (54) 

c(T) may be computed by numerical integration and is plotted in Fig. 4. 

In the intermediate region p ^ 1 F(p, F) must be computed numeri- 
cally. In Fig. 5 we have plotted F(p, F) as a function of p for various 
values of T and in Fig. 6 we have the complete spectra R(a) for the 
same values ofF. 
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FIGURE CAPTIONS 


Figure 1 
Figure 2 


Figure 3 


Figure 4 
Figure 5 
Figure 6 


Angles involved in the scattering process as viewed in the 
electron rest frame. 

Integration boundaries of the variable £ drawn as a function 
of p ' a/a l . (For illustration only - not an accurate plot.) 

a. e <p/[ 1 + y (1 + /3)] 

b. iij [l + y ( 1 + p ) j <e</3 
C. (3 < £ 


. Comparison of approximate and exact expressions for scat- 
tered spectrum from mono-energetic electrons. Initial 
photon energy a x = 10’ 6 . Electron energy is y and D 
= Max [approx/exact - l] 


a. y - Z, 

D 

= .54 

b. y = 9, 

D 

= .024 

00 

rH 

II 

• 

U 

D 

= .0055 


Plot of C(T) versus r. 

. Plot of F(p, T) as a function of p for T = 2, 2.5, and 3. 
. Plot of R(a) as a function of a for T = 2, 2.5 and 3. 
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REGION 1 





REGION 1 


N (ARBITRARY UNI 



Figure 3a ^ 




Figure 3c 
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